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ABSTRACT 

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the 
Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 tt imaging data and 
confirmed with follow-up wide-held photometry from the Large Binocular Cameras. The stellar sys¬ 
tem, with an absolute magnitude of My = —1.8 ± 0.5, a heliocentric distance of SO/l^kpc, and a 
half-mass radius of 34lg pc, shows remarkable similarity to faint, nearby, small satellites such as Will- 
man I, Segue 1, Segue 2, and Bootes 11. The discovery of Laevens 2/Triangulum II further populates 
the region of parameter space for which the boundary between dwarf galaxies and globular clusters 
becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, 
whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar 
structures. 

Subject headings: Local Group — Milky Way, satellites, streams: individual: Laevens 2/Triangulum II 


1. INTRODUCTION 

The last couple of decades saw the discovery of nu¬ 
merous satellites in the Milky Way ( MW) halo. Whi le 
the Sloan Digital Sky Survey (SDSS (jYork et al.ll20nfl[ )) 
satellite discoveries have provided us with greater ob¬ 
servational constraints in our backyard, especially to 
understand the faint end of galaxy format ion in the 
prefe rred cosmological paradigm of ACDM (|Belokurovl 
I20I3D . they have also led to debates about the na¬ 
ture of the faintest satellites (jCilmore et al.l l2007tl . It 
has become apparent that the previously clear distinc¬ 
tion between the compact globular clusters (GCs) and 
the brighter, more extended, and dark-matter domi- 
nated dwarf galaxies (DC s), blurs out for faint systems 
(jWillman fc Stradeil 1201^ . T his is exemplified by the 
discoveries of W illman I (WilI:IWillman et al.ll2005fl and 
Segue 1 (Segl; iBelokurov et ah '2007D. fo llowed up by 
those of Bootes II (BooI I; Walsh et al.ll200^ . and Segue 2 
(Seg2: IBelokurov et al.ll2009D . all nearby satellites within 
25-45 kpc, and just slightly larger than extended outer 
halo GCs. At the same time, these systems are fainter 
than most GCs and all the other DGs. Theoretical expec- 
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tations show that these objects could well be the faintest 
DGs and that tens or hundreds of DGs w ith these prop¬ 
erties could populate the Milky Way halo (|Tollerud et al.l 
120081 : iHargis et al.ll2014l) . As of yet, just two objects have 
been found in the PSl survey (jLaevens et al.ll2?)T^ . re- 
inforcing the tensio n between theory and observations 
(jKlvDin et al.lllOO^ .Onlv the closest DG s would be de¬ 
tected with current pho tometric surveys (|Koposov et al.l 
120071 : IWalsh et al1l2009D . Spectroscopic studies do show 
that the faint systems found so far are dynamically hotter 
than their mere stella r content would impR, hinting that 
they are indeed DGs (iMar tin et al.||2007[ IWillman et al.l 


l2011l : [simon et al1l20III : lKirbv et al.ll20I,3jl . However, the 

low velocity dispersion of these satellites (< 4kms“^), 
combined with the possibl y large impact of binaries 
(jMcConnachie fc C6t^20T(i^ . the complexity of disentan¬ 
gling member stars from foreground contaminants, and 
the overall dimness of their member stars renders any 
definite conclusion difficult. 

Here, we present the discovery of another faint MW 
satellite, Laevens 2/Triangulum 110 , with very similar 
photometric properties to Will, Segl, BooH, and Seg2. 
The new system was found in our ongoing effort to mine 
the Pan-STARRS 1 (PSl) 37r survey for localized stel¬ 
lar overdensities. This letter is structured as follows: in 
section 2, we describe the PSl survey along with the de¬ 
tection method that led to the discovery. We continue 
by discussing follow-up imaging obtained with the Large 
Binocular Cameras in section 3. We discuss the nature 
of the satellite and its implication in section 4. In the 
final section, we summarize and conclude our results. 

In this paper, ma gnitudes are dereddened using the 
iSchlegel et al.l (|1998ll maps, adopting the extinction co- 


® In the absence of spectroscopic confirmation, we wish to remain 
agnostic about the nature of this object and therefore propose a 
double name. For future reference in this paper, we abbreviate to 
Lae 2/Tri 11. 
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Figure 1. : Left: The combined PSl-LBC CMD of all sources within the central 2r^ region of Lae 2/Tri 11. The single epoch PSl 

photometry was used at the bright end (ipi q < 19, squares), with LBC photometry supplementing the faint end 2pi,o > 19 (large dots). 
The orange dashed line indicates the separation between the LBC and PSl data. The red box highlights the clear main sequence of the 
stellar system, the blue box indicates two possible HB stars and the green box identifies likely blue stragglers. Right: Spatial distribution 
of all sources corresponding to the CMD on the left. Large dots correspond to the stars falling within the red CMD box in the left panel 
and show a clear overdensity. The two blue stars indicate the possible HB stars whereas the red ellipse corresponds to the region within 
the favored two half-mass radius of the system, as inferred below. 


efficients of iSchlaflv &: Finkbeined (|2011l) . A heliocentric 
distance of 8 kpc to the Galactic center is assumed. 

2. THE 37r PSl SURVEY AND DISCOVERY 

With a spatial extent encompassing three quarters 
of the sky {S > —30°), PSl (K. Chambers et ah, in 
preparation) gives us an unprecedented panoptic view 
of the MW and its surroundings. Over the course of 
3.5 years, the 1.8m telescope, equipped with its 1.4- 
gigapixel camera covering a 3.3-degree field of view, has 
collected up to four exposures per year in each of 5 bands 
(ffpii’piipi-Zpiypi; iTonrv et al.l 1201^ . A photometric 
catalogue is autom atically generated with the Image 
Processing Pipeline (lMagnieHl2006l . I2007t iMagnier et al.l 
UMl), once the individual frames have been downloaded 
from the summit. The preliminary stacked photometry 
used in this paper has a gpi depth (23.0) that is com¬ 
parable to SDSS y-band depth and rpi/ipi observations 
that reach ~ 0.5/^ 1.0 magnitude fainter : 22.8, 22.5 for 
r and i respectively (jMetcalfe et al.l[20T^ . 

Inspired by past searches for small stellar overdensities 
in MW and M31 surveys, we apply a convolution tech¬ 
nique (Laevens et ah, in preparation), s uccessfully used 
to find new GCs and D Gs in the SDSS (|Koposov et al.l 
l2007t IWalsh et alll2009l l. In a nutshell, we build a mask 
in (r — i,i) color-magnitude space to isolate potential 
metal-poor, old, and blue member stars that could be¬ 
long to a MW satellite at a chosen distance. This mask 
is applied to star-like sources in the stacked PSl photo¬ 
metric catalog. We then convolve the distribution of iso¬ 
lated sources with two Gaussian spatial filters: a positive 


Gaussian tailored to the size of the overdensities we are 
searching for (2', 4', or 8') and a negative Gaussian with 
a much larger kernel (14', 28', or 56'), to account for the 
slowly-varying contamination of sources that fall within 
the color-magnitude mask. By convolving the data with 
the sum of these two (positive and negative) filters and 
accounting for the survey’s spatial incompleteness on the 
arcminute scale, we obtain maps tracking stellar over- 
and under-densities in PSl. We convert these density 
maps into maps of statistical significance by comparison 
with the neighboring regions after cycling through dis¬ 
tances and filter sizes. This procedure already l ed to 
the discovery of Laevens iIB (ILaevens et al.l 120141) also 
discovered concomitantly a s Grat er within the ATLAS 
survey by iBelokurov et al.l (j2014D . The new satellite, 
Lae 2/Tri 11, is located ^ 20° East of M31 and appears 
as a 5.2tT detection, only slightly higher than our sig¬ 
nificance criteria of Soru tailored to weed out spurious 
detections. 

3. FOLLOW-UP 

To confirm the nature and the properties of 

Lae 2/Tri 11, follow-up imaging was obtained with the 
Large Binocular Gamera (LBG) on the Large Binocular 
Telecope (LBT), located on Mount Graham, USA dur¬ 
ing the night of October 17-18 2014. With its 23' x 25' 

Following th e naming convention established in 

IBianchini et al] 11201511 

These also include a check that potential detections do not 
also correspond to a s ignificant overdensity of background galaxies 
UKodosov et al]l2nn7ll . 
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Figure 2. Left: The CMD of Lae 2/Tri II within 2r^, with best-fit isochrones overplotted. Both isochrones have an age of 13 Gyr and 
are shifted to a distance modulus of 17.3. The red/blue isochrones have metallicities [Fe/H] = —2.19 a nd [Fe/H] = —1.80 , respectively. 
Right: The same CMD of Lae 2/ Tri II with fou r metal-poor fiducial isochrones from GCs observed in PSl UBernard et al.ll2014l) . dereddened 
assuming the reddening values of (IHarriall201CII) and shifted to match the observed features. The green fiducial (NGC 7078; [Fe/H] = —2.37), 
shifted to a distance-modulus of 17.5 best represents the features of the new stellar system, with the other fiducials i.e. red: NGC 4590 
([Fe/H] = —2.23), blue: NGC 6341 ([Fe/H] = —2.31), orange: NGC 7099 ([Fe/H] = —2.27) appearing too red or to blue to accurately 
reproduce the MS and MSTO. NGC 4590 is shifted to a distance-modulus of 17.5, whereas NGC 6341 and 7099 are at 17.3. 



field of view and equipped with 4 CCDs, the LBC are 
ideal to follow-up MW satellites that usually span a few 
arcminutes on the sky. Imaging was conducted in the g 
and i bands, making use of the time-saving dual (binoc¬ 
ular) mode using the red and blue eye simultaneously. 
Six dithered 200s sub-exposures were acquired in each 
band, with a seeing of 1". The field was centered on the 
location of Lae 2/Tri II. 

The images were processed and the photometry 
performed using a version of the CASU pipeline 
(jirwin fc LewisI 120011 ) updated to work on LBC data. 
The instrumental magnitudes were calibrated onto the 
PSl system {gpi and ipi), by compar ison with the PSl 
single epoch data (iSchlaflv et ahllMl^ to derive the rele¬ 
vant color equations. The final LBC photometry reaches 
more than 2 magnitudes deeper than the stacked PSl 
data. The left-hand panel of Figure [T] shows the com¬ 
bined PSl/LBC color-magnitude diagram (CMD) of all 
stars within 2 half-light radii (±2r/j), for which all sources 
brighter than fpi.o = 19.0 are taken from the PSl sin¬ 
gle epoch photometry so as to extend the CMD beyond 
the saturation limit of the LBC photometry. The main 
sequence (MS) of an old and metal-poor stellar system 
is readily visible, with a clear turn off at fpi^o ^ 20.8. 
Given the low density of the MS, the red giant branch 
(RGB) of the stellar system is very likely sparsely pop¬ 
ulated and hidden within the foreground contamination 
of bright stars. However, two blue stars are consistent 
with being blue horizontal branch (HB) stars and are 
highlighted in the CMD. The other panel of Figure [T] 
presents the spatial distribution of sources in the four 


chips, with a red ellipse indicating the satellite’s two 
half-light radii extent (as determined by the structural 
parameter analysis, see section 4). Stars with colors and 
magnitudes consistent with the MS are shown as large 
dots and reveal a clear spatial overdensity. Likely blue 
straggler stars are identified and given by the green tri¬ 
angles. The two blue stars correspond to the potential 
HB stars; their location close to the center of the stellar 
overdensity supports them being member HB stars. 

4. PROPERTIES OF THE STELLAR SYSTEM 

Since the RGB and HB are so sparsely populated, 
an investigation into the presence of member RR Lyrae 
stars in the multi-epoch PSl data unsurprisingly led to 
no candidate from which to derive a distance estimate. 
However, due to the well defined MS and MS turnoff 
(MSTO) at ipi,o 20.8, a reliable distance estimate can 
nevertheless be determined through a comparison with 
isochrones and fiducials by eye (Figure [2]). We assess the 
stellar system’s metallicity, age, and distan ce modulus by 
first cycling through Parsec isochrones (jBressan et al.l 
I2012D for a metallicity range —2.2 < [Fe/H] < —1.3 
{Z = 0.0001 to 0.0007, assuming Zq = 0.152) and an 
age between 9 and 13 Gyr. We investigate the effects that 
various metallicities and ages have on the distance deter¬ 
mination, by cycling through distance-modulus steps of 
0.1 between 16.9 and 17.5. The isochrone which best rep¬ 
resents the CMD features is a metal-poor, old isochrone 
([Fe/H] = —2.19 and age of 13 Gyrs), for a distance mod¬ 
ulus of 17.3. We further strengthen these conclusions by 
comparing the CMD of Lae 2/Tri H with the fiducials 
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Figure 3. Left: Probability distribution functions for the ellipticity (e), the position angle (0), the angular and the physical half-mass 
radii (r^) of Lae 2/Tri II (from top left to bottom right). Middle: Comparison between the favored radial distribution profile (full line) 
and the data, binned according to the preferred structural parameters (dots). The error bars assume Poissonian uncertainties and the 
dashed line represents the field density. Right: Probability distribution function of the absolute magnitude of Lae 2/Tri II in the V band 
{My = -1.8 it 0.5). 


from 13 GCs and 3 Open Clusters of varying metallicity , 
derived directly from the PSl data (iBernard et al.l[20l3) . 
The 4 most metal-poor GCs of the sample provide a good 
fit to the MS and MSTO of Lae 2/Tri II provided they 
are shifted to distance moduli in the 17.3-17.5 range. 

Combining these two analyses by averaging the six 
best-fits fiducials and isochrones, we therefore conclude 
that Lae 2/Tri II is old and metal-poor, and is located 
at a distance modulus of 17.4to,i, which translates to a 
heliocentric distance of 30^2 kpc, or 36^2 kpc from the 
Galactic center. In both cases (isochrone and fiducial), 
the HB of the favored track also overlaps almost exactly 
with the two potential HB stars. Please note that the 
uncertainty in the distance measurement is propagated 
through for the derivation of the structural parameters, 
further detailed in the next few paragraphs. We also 
draw to the reader’s attention that a more involved anal¬ 
ysis of ‘CMD-fitting’ would likely enhance the quality of 
the distance measurement; however, the limitation of the 
field of view prevent us from obtaining a large enough 
sample of background stars to perform such an analysis. 

We derive the structural parameters of Lae 2/Tri II 
by using a mod i fied v ersion of the technique described 
in iMartin et al.l (1200811 . updated in Martin et al. (2015, 
submitted). The updated technique allows for a Markov 
Chain Monte Carlo (MCMC) approach, based on the 
likelihood of a family of exponential radial density pro¬ 
files (allowing for flattening and a constant contamina¬ 
tion over the field) to reproduce the distribution of the 
system’s MS stars. The parameters of the model are: the 
centroid of the system, its ellipticitj{3, its position an¬ 
gle (defined as the angle of the major axis from North to 
East), its half-mass radiusQ, fh, and the number of stars, 

The ellipticity is here defined as I — b/a with a and b the 
major and minor axis scale lengths, respectively. 

Note that, assuming no mass segragation in the system, the 


Table 1 

Properties of Lae 2/Tri II 


a (J2000) 

02:13:17.4 

<5 (J2000) 

-1-36:10:42.4 

£ 

141.4° 

b 

-23.4° 

Distance Modulus 

d d 
+ 1 

Heliocentric Distance 

30/2 kpc 

Galactocentric Distance 

36/2 kpc 

Mv 

-1.8±0.5 

[Fe/H] 

~ -2.2 

Age 

^ 13 Gyr 

E{B - VY 

0.081 

Ellipticity 

0.2l/°;i( 

Position angle (from N to E) 

56+is° 

rh 

0 Q-l-1.1/ 

'^•^- 0.9 

rh 

34t« pc 


“ from ^chlege^^t_^ (T998) 


N*, within the chosen CMD selection box. Although it 
is located in the Triangulum constellation, Lae 2/Tri II 
is so far from both M31 (~20 degrees) and M33 (^10 
degrees) that any contamination by M31 or M33 stellar 
populations is vanishingly small and does not impact our 
results. The resulting probability distribution functions 
are presented in the left-most panels of Figure |3] for the 
most important parameters and summarized in Table 1. 
Figure [3] also shows a favorable comparison of the pre¬ 
ferred exponential profile with the data binned following 
the preferred structural parameters. 

It should however be noted that the properties of 
Lae 2/Tri II as observed by the LBC could be slightly 
biase d by the low contras t of the stellar overdensity. In¬ 
deed, iMuhoz et al.l (1201211 show that satellite properties 

half-mass radius is equivalent to the half-light radius. 
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Figure 4. The distribution of MW satellites in size—magnitude space. GCs are shown as squares, DGs are shown as circles and Lae 2/Tri II 
is represented by the large star symbol. The color scale indicates the ellipticity of the various satellites. Lae 2/Tri IPs ellipticity and half¬ 
mass radius show very similar values to those of the four satellites: Segl, Seg2, BooII and Will. Finally, w e also indicat e the recently 
disco vered MW sat e llite, Laevens 1/Crater given by the large triangle. The data for the GCs were taken from IHarrisI II2010I 1 and the DGs 
from IMcConnaclilS 112012 ). 


are most accurately measured when the central density 
of stars relative to that of the background is larger than 
20, which is not the case here. Deeper data would be 
necessary to strengthen our size measurement. 

To determine the absolute magnitude of the stel¬ 
lar system, we foll ow the same proced ure we used for 
Laevens 1/ Crater (iLaevens et al.ll2014D as was initially 
described in iMartin et al.l (j2008D . After drawing a value 
of N* from the structural parameter cha i n, we sample the 
CMD of the best-fitting iBressan et al.l (j2012D isochrone 
(see Figure 12), with its associated luminosity function 
and photometric uncertainties, until it contains N* stars 
in the CMD selection box used for the structural param¬ 
eter analysis. Adding up the flux of all the stars drawn 
in this artificial CMD yields the absolute gpi- and ipi- 
band magnitudes of Lae 2/Tri II {Mg = —1.7 ± 0.5 and 
Mi = —2.1 ± 0.5), which converts to My = —1.8 ± 0.5. 
This technique has the benefit of accounting for the effect 
of sampling such a small population of stars may have on 
the determi nation of the system ’s magnitude (i.e. CMD 
‘shot-noise’; IMartin et al.H200^ . 

5. DISCUSSION AND CONCLUSION 

We have presented the discovery of a new MW satel¬ 
lite, Lae 2/Tri II, discovered within the PSl Stt data 
and confirmed from deep and wide LBC follow-up. Lo¬ 
cated at a heliocentric distance of SOt^kpc, this system 
is very faint {My = —1.8 ±0.5), old 13 Gyrs), metal- 
poor ([Fe/H] ~ —2.2), small (34+1 pc), and mildly el¬ 
liptical (0.21+g'2[). Figure 0] places Lae 2/Tri II in re¬ 
lation to other MW GCs and DGs. This new system’s 
magnitude and half-mass radius are very similar to the 
properties of the faint satellites Segl, Seg2, Will, and 


BooII, which were all recently discovered in the SDSS. 
Ultimately, high quality spectroscopic follow-up and an 
assessment of its dynamics are necessary to confirm the 
nature of this new satellite. However, its similarity in dis¬ 
tance, size, absolute magnitude, age, and metallicity to 
those of Will, Segl, BooII, and Seg2, that all have larger 
velocity dispersion than implied by their t iny stellar mass 
(IMartin et al.l [20 ^ iSimon et all 120111 : iWillman et al.l 
120111: iKirbv et al.l 120131 ) hints that Lae 2/Tri II could 
well be another one of these systems that appear to pop¬ 
ulate the faint end of the galaxy realm. 

It is also worth noting that the location of Lae 2/Tri II, 
{e,b) = (141.4^-23.4°), - 20° East of M31, places it 
within t he Triangulum-Andromeda stellar structure(s) 
(TriAnd: iMaiewski et al.ll2004l: iR.ocha-Pinto et ai1l2004t 

iSheffield et al.l 120141) . Although this MW halo stellar 
overdensity is very complex, with e vidence for multiple 
substructures (iBonaca et al.l 120121 : IMartin et al.l 120131 : 
IMartin et al.l I2014D . it~ spans a large enough distance 
range to encompass Lae 2/Tri II (~ 15 — 35kpc). A 
recent spectr o scopic study of stars within TriAnd by 
iDeason et al.l (120141) confirmed t hat, as initially pro- 
posed bv iBelokurov et all (I2009D . Seg2 is also likely 
embedded within it and follows a systematic trend of 
these faint satellites being part of MW halo stellar 
streams. Segl has been proposed t o be tied to the Or¬ 
phan Stream (jNewberg et al.l[2010D . though differences 
in abun dance patterns between both have also been ob¬ 
served (I Vargas et al.l 120131 : iCasev et al.l I2014D . Simi¬ 
larly Boon’s distance and radial velocity are compatible 
with it being part of the Sagittarius stream (|Koch et al.l 
[2?)00l) . whereas high resolution abundance measurements 
for BooII stars question this association (jKoch fc Richl 
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I2014f) . It remains possible, however, that the small stel¬ 
lar systems were satellites of the larger, now disrupted 
progenitor of these stream, thereby alleviating the need 
for them to share similar abundances. In this context, it 
is particularly interesting that Lae 2/Tri II is situated on 
the linear extrapolation of the Pan-Andro meda Archae¬ 
ologi cal Survey (PAndAS) MW stream (|Martin et al.l 
1201^ . 10° beyond the PAndAS footprint where this 
dwarf galaxy remnant was discovered. The stream and 
satellite are not aligned; however, the uncertainties on 
the position angle and ellipticity are not conclusive in 
ruling this out. Here as well, spectroscopy is necessary 
to derive the systemic velocity of Lae 2/Tri II and confirm 
it is compatible with the global motion of Triangulum- 
Andromeda and, in particular, with the velocity of the 
PAndAS MW stream. Follow-up will help reinforce or 
disprove such a hypothesis. 

Although PSl is only slightly deeper than the 
SDSS, the extra coverage provided by its Stt footprint 
leaves hope for more discoveries of faint objects like 
Lae 2/Tri 11. Building up the statistics of these systems 
through more discoveries in current (PSl) a nd (DES; 
iThe Dark Energ y Survey C ollaborationI 120051) or future 
surveys ILSST: lTvsonll2002h is essential if we are to un¬ 
derstand the true nature of these incredibly faint stellar 
systems that can only be found in the MW surroundings. 
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